Anticipated synchronization in coupled 
chaotic maps with delays 



Cristina MasoUer^, Damian H. Zanette^ 

^Instituto de Fisica, Facultad de Ciencias, Universidad de la Repuhlica, Igud J^225, 

Montevideo 114-00, Uruguay 

^Consejo Nacional de Investigaciones Cientificas y Tecnicas, Centro Atomico 
Bariloche and Instituto Balseiro, 84OO Bariloche, Rio Negro, Argentina 



Abstract 

We study the synchronization of two chaotic maps with unidirectional (master-slave) 
coupling. Both maps have an intrinsic delay ni, and coupling acts with a delay n2- 
Depending on the sign of the difference ni — n2, the slave map can synchronize 
to a future or a past state of the master system. The stability properties of the 
synchronized state are studied analytically, and we find that they are independent 
of the coupling delay n2- These results are compared with numerical simulations of a 
delayed map that arises from discretization of the Ikeda delay-differential equation. 
We show that the critical value of the coupling strength above which synchronization 
is stable becomes independent of the delay ni for large delays. 
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1 Introduction 



Time-delay systems liave attracted a lot of attention in recent years, in part 
due to the fact that multistability, i.e. the coexistence of multiple attractors, 
is a common occurrence when the delays are large -typically, much larger 
than the response time of the system [1]. Interest in multistability arises be- 
cause multistable systems play a key role in pattern recognition processes [2] 
and memory storage devices. By choosing appropriate initial conditions, pre- 
scribed periodic solutions can be stored as oscillatory patterns of a time-delay 
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system [3-5] . Synchronization of chaotic time-delay systems has also received 
attention, since it has potential applications to secure communications [6,7]. 
Perez and Cerdeira [8] have shown that, in low-dimensional chaotic systems, 
a hidden message can be unmasked by the dynamical reconstruction of the 
chaotic signal using nonlinear dynamical methods. Encrypting a message in 
the chaotic output of a time-delay system has the advantage that the dynam- 
ics is in this case high-dimensional (the dimension increases linearly with the 
delay [9]) but, in spite of this fact, synchronization can be achieved by trans- 
mitting a single scalar signal. Unfortunately, this method is not as secure as 
initially expected, since it has been shown that by using a special embedding 
space, the delay time can be identified, and the message can be successfully 
unmasked [10]. Coupled oscillators with time delays in the coupling, which 
represent interactions being transmitted at finite speed, have been extensively 
studied as well [11-13]. 

Recently, a new effect of delayed coupling was reported by Voss [14,15], who 
showed the existence of an anticipating synchronization regime. In this regime, 
the slave system becomes synchronized to the chaotic future state of the master 
system. Anticipation occurs when the coupling is delayed, and results from 
the interplay of memory effects and relaxation mechanisms. Numerically, this 
regime was found in coupled semiconductor lasers with optical feedback [16]. 

In this paper we investigate the existence and stability of anticipated synchro- 
nization in coupled time-delay maps, which have the advantage of allowing 
for analytical calculations. In the next section, we introduce a system of two 
time-delay coupled maps in a master-slave configuration, and present ana- 
lytical results on the stability of anticipated and retarded synchronization for 
generic maps. In Section 3, we apply the results to a delay map that arises from 
the discretization of the Ikeda delay-differential equation. Section 4 presents 
a summary and the conclusions. 



2 Master-slave coupled delay maps 

We consider a one-dimensional map of the form 

Xn+i = hxn + /(a;„_nj. (1) 

For \b\ < 1, the first term in the right-hand side represents a relaxation mech- 
anism. Under its sole action Xn would asymptotically vanish. This relaxation, 
however, competes with the effect of the nonlinear function f{x), which has 
the form of a time-delayed feedback with delay rii. 
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The (master) map (1) is used to partially drive the evolution of a new (slave) 

system The dynamics of this system is, in principle, the same as for Xn, 
except that a part of the nonlinear component is replaced by the evolution of 
Xn, namely 

Un+i = &yn + (1 - r])f{y n—ni n—n2 )■ (2) 



The parameter vj G [0, 1] measures the strength of the driving action. Note that 
Xn enters the dynamics of Un with a delay ^2- This form of unidirectional time- 
delay coupling is not only an extension to maps of the configuration proposed 
by Voss for continuous-time systems [10], but also generalizes the master-slave 
interaction by introducing a parameter (rj) that controls its strength. Voss's 
configuration corresponds to the extreme value t] = 1. 

The structure of Eqs. (1) and (2) suggests that, if driving is effective enough, 
Hn may become synchronized to the driving signal, in the form 

Un •^n—n2+n\ • (3) 



If this synchronization is realized, the dynamics of y„ coincides with that of 
Xn up to a time shift rii — n2- If ni > n2 the slave system anticipates the 
evolution of the master system, and we have anticipated synchronization. On 
the other hand, if rii < n2 the synchronization is retarded. 

Linear stability of the synchronized state (3) can be analyzed in the usual way. 
Taking yn = Xn-n2+ni + Syn, the evolution of the deviation dyn in the limit of 
infinitesimal deviations is readily derived from (2): 

5yn+l = bSyn + (1 - ri)f'{Xn-n2)Syn-ni: (4) 



with /' = df/dx. Note that, if \b\ < 1 and for t] = 1, the deviation decreases 
exponentially with time and synchronization is therefore stable. In view of 
well-known properties of synchronization in coupled maps without delays, it 
is possible to advance that, if the dynamics of x„ is chaotic, the synchronized 
state (3) will be stable above a certain threshold rjc < 1. On the other hand, 
if Xn is nonchaotic, synchronization should be stable for any rj > 0. 

For Ui 7^ 0, Eq. (4) can be formally integrated by introducing a linear (ni-|-l)- 
dimensional map for a variable z„ = {z^, z^^, . . . , z^^^) , with z^ — Syn-k- This 
equivalent map is given by 

Zn+l = MnZn, (5) 
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where the elements of the matrix M„ read 



b 



if i 



J = 0, 



(1 - V)f{x, 



) if i 



(6) 



1 



Hi 



j + 1 (i = 0,...,ni-l), 







otherwise. 



The deviations given by Eq. (4) can then be obtained from the solution of the 
equivalent map, 



Thus, synchronization is linearly stable if all the eigenvalues of the evolution 
matrix [/„ vanish for n — > cxd. For rj = 1 the eigenvalues are (6", 0, . . . , 0). For 
Tj < 1, on the other hand, Un and its eigenvalues must be calculated numer- 
ically, since they involve the successive states of the variable x„, cf. Eq. (6). 
It is interesting to point out that stability is determined by the asymptotic 
properties of f/„ in the limit n — > oo, where the driving system a;„ has exhaus- 
tively explored the accessible region of phase space. In this limit, thus, the 
delay n2 of Eq. (2) plays no role in the determination the eigenvalues of Un- 
In consequence, n2 is irrelevant to the stability of the synchronized state. 



3 Application to the Ikeda delay map 

In this section we present numerical results corresponding to maps (1) and (2), 
with = asinx. With this choice of f{x), we refer to Eq. (1) as the Ikcda 
delay map, since it arises from the discretization of the Ikeda delay-differential 
equation [1,17], 

X — —X + asmx{t — t). (8) 

This equation was introduced to describe the dynamics of an optical bistable 
resonator, where the finiteness of the speed of light is relevant, and therefore 
require a description that incorporates explicitly the round-trip time of light 
in an optical cavity. Physically, x is the phase lag of the electric field across 
the resonator, a is the laser intensity injected into the system, and r is the 
round trip time of the light in the resonator. 

Writing x = [x{t + At) — x{t)]/ At, t = nAt, and r = niAt, the Ikeda delay 
map is obtained from Eq. (8) with b = 1 — At, a = a At, and Xn = x{t). 



z 



'n 



UnZo = Mn-lMn-2 ■ ■ ■ MiMqZq. 



(7) 
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Fig. 1. Time series of the maps (1) solid line) and (2) dashed line) with 
f{x) = asinx, for a = 3, 6 = 0.9, r) = 0.833, ni = 15, and n2 = 0. The dotted line 
stands for the difference Xn+m-n^ ~ Vn- 

The Ikeda delay map is not to be confused with the discrete Ikeda map [1,18], 
which is obtained from (8) by discretizing time in units of r in the singular 
limit where the delay-to-response time ratio diverges. 

In the following, we focus the attention on the case b = 0.9. Our numerical 
calculations are restricted to random initial condition in (— tt, tt) both for the 
master and the slave system. In Fig. 1 we illustrate anticipated synchroniza- 
tion, for ni = 15 and n2 = 0. The slave system (dashed line) anticipates in 
15 steps the state of the master system (solid line). After a transient, the dif- 
ference Xn+ni-n2 ~ Vn (dotted line) decays to zero. We have verified that this 
behavior is independent of the value of n2- 

The critical value rjc of the coupling strength, above which the synchronized 
state is stable, has been calculated from the analytical results of the previous 
section by means of a numerical evaluation of the matrix Un in Eq. (7). Since 
the delay n2 is irrelevant to this calculation, we can take ^2 = in Eq. (6). 
Prom a given initial condition for the master system, the matrix C/„ is evaluated 
for successive values of n. In order to avoid the calculation of its eigenvalues, 
which is specially troublesome for large rii, Un is multiplied at each step by 
a randomly generated vector, u^, of unitary modulus. If, for a given value of 
77, the modulus of the product remains above a certain upper threshold H for 
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Fig. 2. Critical value of the coupling strength above which the synchronized state is 
stable, calculated by means of a numerical evaluation of the matrix Un in Eq. (7), 
as explained in the text. 

sufficiently large n, |t/nUn| > H iov n > nmax, synchronization is considered 
unstable and r]c > rj. If, on the other hand, |C/riU„| < h ior n > rimax, where h 
is a suitable lower threshold, synchronization is considered stable and rj^ < V- 
With this criterion, rjc can be found by decimation within the interval (0, 1) 
up to a certain previously fixed precision. In our calculations, we have taken 
H = IQio, h = 10-^°, and n^nax ~ 10^ 

Figure 2 shows function of the delay ni and for several values of the 

parameter a. Note that for small rii -specifically, for rii = 1- the critical value 
r]c can vanish, which indicates that the master evolution is nonchaotic. Notice, 
furthermore, that for some values of a and rii — 1 we have plotted more 
than one value of rjc- This is due to the fact that, in such cases, the master 
system is multistable -with, typically, two chaotic attractors. The value of 
rjc depends thus on the initial conditions for both Xn and y„. Whereas the 
behavior is quite irregular for small ni, for larger delays the critical coupling 
becomes practically independent of rii. Though it has not been possible to 
prove this analytically, we conjecture that rjc approaches a finite limit below 
unity as ni oo. The dependence on a in the investigated interval is also 
moderately weak. As expected, since the dynamics of the Ikeda delay map 
becomes more irregular as a grows -i.e. the Lyapunov exponent is higher- rjc 
increases accordingly. 
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Fig. 3. Similarity function, Eq. (9), for the same parameters as in Fig. 1 
(ni - n2 = 15). 

The degree of anticipated of retarded synchronization can be quantified by 
calculating the similarity function Sj, defined as [19]: 



If Xn and Un are independent time series with similar mean value and disper- 
sion, the average of their square difference is of order 2{xD ?a 2(|/^), and thus 
Sj ^ -\/2 1.4. If, on the other hand, there is perfect anticipated or retarded 
synchronization, the difference Xn+ni-n2 ~yn vanishes, and Sni-n2 = 0- Fig- 
3, the similarity function is shown for the same parameters as in Fig. 1. 

Figure 4 shows min(5') -i.e. the minimum value of Sni-n2 recorded during 
sufficiently long realizations of the evolution for a given set of parameters- in 
the (ni, 7/)-plane, with all the other parameters as in Fig. 1. The region where 
anticipated synchronization occurs, min(S') ~ 0, is clearly visible for large 
values of rj. The position of its boundary is in good qualitative agreement 
with the values of rjc shown in Fig. 2 for a = 3. 
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Fig. 4. Minimum of the similarity function in the parameter space (77, ni), with all 
the other parameters as in Fig. 1. The minimum was calculated by averaging over 
50 trajectories with different initial conditions. 

4 Conclusion 

We have extended the study of anticipated synchronization, advanced by Voss 
for unidirectionally coupled differential equations with time delays [14,15], to 
delayed coupled chaotic maps. While the nature of anticipated synchronization 
of maps and differential equations is the same, delay discrete-time dynamics 
admits an analytical treatment which cannot be carried out for continuous- 
time systems. In fact, ordinary differential equations with finite time delays 
constitute an infinite-dimensional problem. On the other hand, since times 
delays in maps must be discrete, the dimensionality of the problem remains 
finite. 

Taking advantage of this situation, we have analytically studied the stability 
of anticipated and retarded synchronization in a generic master-slave config- 
uration. In the absence of coupling, master and slave dynamics are identical 
and involve an intrinsic delay rii. Coupling consists in the replacement of a 
part of the slave dynamics by that of the master system, with a delay 77.2. We 
have shown that the stability of synchronization is independent of n2. The 
structure of the linearized problem, Eqs. (4-7), suggest meanwhile a strong 
-though not transparent- dependence on the Lyapunov exponent of the mas- 
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ter system, as expected. In practice, the linearized problem has to be treated 
numerically, but it only involves the realization of the master system and the 
successive application of the linear-evolution operator, being thus a purely 
algebraic process. 

These results have been applied to the Ikeda delay map, which derives from 
the application of the Euler integration scheme to the Ikeda delay-differential 
equation. We have calculated the critical coupling intensity above which syn- 
chronization is stable, as a function of the delay ni and of a parameter a 
that controls the chaotic dynamics of the map. It is found that, whereas for 
small values of rii the critical coupling can vary considerably -due to the ir- 
regular appearance and disappearance of chaotic and nonchaotic attractors- 
the dependence for large rii is much smoother. In fact, as rii grows, the crit- 
ical coupling seems to approach a constant value. The dependence with the 
dynamical parameter a is also moderate in the considered range. 

In this work we have focused the attention on exact anticipated synchroniza- 
tion. However, it has been previously shown that approximate anticipated syn- 
chronization is possible in coupled differential equations, even in the absence 
of intrinsic delays [15]. The study of approximate anticipated synchronization 
in coupled maps constitutes therefore a line for future work. In particular, it 
would be interesting to investigate in detail the connection between the degree 
of synchronization, and the irregularity of the dynamics -as measured by the 
Lyapunov exponents. For delay differential equations, the number of positive 
Lyapunov exponents and the fractal dimension increase linearly with the de- 
lay, while the metric entropy remains roughly constant [9,20,21]. We therefore 
conjecture that the metric entropy might be a good indicator of the possibility 
of synchronizing with anticipation, and thus predicting, chaotic dynamics. 
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